Fidelity of photon propagation in electromagnetically induced transparency in the 

presence of four- wave mixing 
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We study the effects of the four-wave mixing (4WM) in a quantum memory scheme based on 
electromagnetically induced transparency (EIT) . We treat the problem of field propagation on the 
quantum mechanical level, which allows us to calculate the fidelity for propagation for a quantum 
light pulse such as a single photon. While 4WM can be beneficial for classical, all-optical information 
storage, the quantum noise associated with the signal amplification and idler generation is in general 
detrimental for a quantum memory. We identify a range of parameters where 4WM makes a single 
photon quantum memory impossible. 

PACS numbers: 42.50.Gy, 42.50.Ct, 03.67.Hk, 42.65-k 



I. INTRODUCTION 

A reliable quantum memory for photons is one of the 
essential ingredients for quantum networks and optical 
quantum computing. There have been several proposals 
for photon storage, which fall into three main categories: 
gradient echo techniques [1], far detuned Raman systems 
[2], and electromagnetically induced transparency (EIT) 
[3]. In all these schemes the storage of single photons 
occurs by nrapping the quantum state of photons onto a 
long lived atomic excitation. 

In this paper we will concentrate on the EIT based 
scheme, where one uses a strong control field to couple an 
incoming signal pulse to the atomic spin coherence result- 
ing in the common propagation of both as a dark state 
polariton. By adiabatically switching off the control field 
the signal field is mapped on the spin coherence and later, 
after some storage time, is retrieved by switching on the 
control field. Since its theoretical proposal [4, 5] and the 
first experimental realizations [6, 7] there has been a large 
development of EIT based quantum memories (QM), e.g. 
successful implementation in hot gases [8], in cold gases 
using magneto-optical traps (MOT) [9] or optical lattices 
[10, 11], as well as solid systems such as rare-earth doped 
crystals [12]. E.g., using EIT memory in gases, weak 
coherent pulses have been stored in hot Rubidium gas 
with storage times of T s = lms and storage efficiencies 
of r] m = 45% [13]. While in a cold Rb MOT, EIT has 
been used to store pulses with storage times of T s = lms, 
with?7 m = 78% [14]. 

The storage efficiency of EIT QM is limited by two 
considerations. First the spatial pulse size L p must fit 



entirely inside of the medium L p 



T p v g 



< L. where T v 



is the pulse duration and v g is the group velocity in the 
medium, otherwise some of the pulse will leak out during 
the storage process and be lost. Secondly, the spectral 
width of the pulse Au> p must be well within the EIT 
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transmission window, Aoj p <C ujeit — \Dv g /L, where 
D = L/L^bs denotes the optical depth of the medium, i.e. 
the ratio of medium length L to absorption length L a bs 
in the absence of EIT. Since the spectral width and pulse 
length are inversely proportional Acj p ~ 1/T P , these two 
requirements compete with each other and both of them 
can only be satisfied at large optical depth VD ^> 1 [15]. 

However with high optical depth, non-linear processes 
start to become important. In particular, a four wave 
mixing process (4WM) is possible in many of the imple- 
mentations of EIT QM, where the control field with Rabi 
frequency f2 and appropriate polarization also acts as a 
far-detuned field with Rabi frequency 51' on the signal 
transition spontaneously generating a new 'idler' field. 
This idler field then moves population into the spin state, 
which is then pumped by the control field to the excited 
state as shown in Fig. 1. The medium is still transparent 
to the signal pulse due to EIT, but now the signal pulse 
also experiences some gain from 4WM. 

The role of 4WM in EIT quantum memories is not 
completely clear. In one experiment, 4WM was sug- 
gested as an explanation of why the storage efficiency 
has tended to saturate to values lower than 50% with 
high D in EIT QM experiments [16]. However, there 
have also been experiments [17] that claim 4WM may be 
useful due to better spatial pulse compression and pulse 
gain. [17] also suggests that with the help of the 4WM 
one may achieve multimode storage, storing not only the 
signal mode, but also the idler mode. However this con- 
jecture was disproved in a more recent experiment [18] 
where it was clearly shown that multimode storage is not 
possible in this system, due to no significant slowing of 
an input idler field, allowing it to escape the medium be- 
fore storage. Therefore, we will not consider the case of 
an input idler field. 

EIT with 4WM has the advantage of signal gain which 
could be used to compensate losses in the medium, natu- 
rally improving the storage of classical signal pulses. But 
the goal of a quantum memory is single photon storage, 
where gain can become a liability since it is always accom- 
panied by additional noise. We therefore address the case 




FIG. 1: Level scheme for EIT memory with four- wave mixing 
process. A double A-scheme, with one A being the signal 
field as and strong control field Q giving standard EIT, and 
a second A far-detuned from resonance made up of the same 
control field acting on the \g) — \e) transition and the idler 
field a i generated by 4WM. 



nal photon form the 



\g) transition, il = fi ' s r Ec 



and n' 



AV s £c 



are the Rabi frequencies of the con- 



trol field with dipole moments ^ eg , /Vg, N is the num- 
ber of atoms in the medium, and V is the volume of the 
medium, ff^(z) are continuous atomic ensemble spin- 
flip operators corresponding to the transition from in- 
ternal state \v) to \fi), and as and ai are dimcnsionlcss 
field operators of the signal- and idler fields, which fulfill 
bosonic commutation relations [a(z),a)(z')} = S(z — z'). 

Qs = VegJ^v, 9i = Ve's^J^y arc the coupling con- 
stants for the field operators as and <2j. 

The equations of motion for the atomic operators are 
given by the Heisenberg-Langevin equations: 

d . 



dt 



\H., 



int,Cr]ij — "fijVij + 5ij 2_^ TUUU + F, 



*J> 



(2) 



of a single photon propagating in an EIT medium with 
4WM, by developing a fully quantum model for pulse 
propagation in Sect. II, which we solve in Sect. III. 
We then analyze the amplification noise in the system in 
Sect. IV resulting from the spontaneous generation of 
idler photons coupling to the signal field. In Sect. V we 
use our results to calculate the memory fidelity of EIT 
with 4WM in an otherwise loss- less medium. In Sect. VI 
we consider additional linear losses and discuss the case 
where 4WM gain completely compensates these losses 
in the medium. Finally, as in EIT memories the fields 
are close to atomic resonances we expand our analysis in 
Sect. VII to incorporate noise associated with population 
decay from excited states. 



II. MODEL 

We can model the EIT 4WM system as an ensemble of 
four level atoms, which interact with the strong control 
field E c coupling the \s) — |e) as well as acting as a far- 
detuned field on the \g) — \e') transition and a weak co- 
propagating signal field Es in resonance with the \g)— |e) 
transition, see Fig. 1, a treatment similar to [19]. Due to 
the additional coupling of the control field to the \g)— \e'} 
transition a new idler field will be generated which be- 
cause of energy and momentum conservation, propagates 
then in the same directions as the other two fields. In this 
approach |e') can be considered either as an additional 
level, or as a representation of the excited state to sim- 
plify the coupling of two fields to the same transition, in 
cither case it will eventually be adiabatically eliminated. 
The interaction Hamiltonian in the dipole, rotating wave, 
and slowly varying envelope approximations is given by: 



Hint = 



™f*z{6» e 
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where A is an effective detuning of the control laser from 
the |e') — \g) transition, S is the detuning of the sig- 



where 7^ are the decoherence rates, 77, are the spon- 
taneous emission rates from |Z) to \i) and Fij arc 5- 
correlated Langcvin noise operators. The evolution of 
the fields is governed by the following propagation equa- 
tions 



d d . 

d-t +C d-z) 
d d . 



(3) 

(4) 



We simplify these equations by assuming the signal and 
idler fields remain weak, such that they can be treated 
perturbatively in the atomic equations. This effectively 
fixes all of the population in the ground state. Therefore, 
the strong control field is not significantly depopulated 
and we can assume fi and fi' are constant. We then adi- 
abatically eliminate |e') leaving four coupled equations: 



idtd-ge = (6„ — 6- i"f gc )(7ge - ga s + ifta gs +iF t 



y< : 



(5) 
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t u gs 



(5 s -6 -ij gs )a 



gs 



a' + 



(d t + cd z ) a s = igNa gei 

Of* 

(d t + cd z ) a\ = -igN—a gs , 



n*a ge +iF gs , (6) 
(7) 
(8) 



where 6s = |S!'| 2 /A is the AC-Stark shift and for sim- 
plicity we take gs = gi = g- The same equations can be 
derived directly from a 3-level model, as done by Phillips 
ct al. [18], with the only difference being an additional 
AC-Stark shift in Eq.(5), (8 S - S) -> (25 s - 5) that can 
not be removed by the choice of dctunings. Since we are 
considering the far detuned regime with A 3> 7 gc , this 
frequency shift S s can be neglected as it will be much 
smaller than the EIT transmission window. 



III. PULSE PROPAGATION 

In the case of constant control field we can easily solve 
Eqs.(5-8). For simplicity we take the single photon de- 
tuning to match the AC-Stark shift, 5 = S s = \Q\ 2 /A, 



and set 7 gs = 0. The solutions in terms of the optical depth D, are given in the frequency domain as: 
I 

a s (D,co) = A(D,u) a s (0,cu) + B(D,u) 4(0, -w) + 5a Sl 
a\(D,u) = -B(D,u>) a s (0,u>) + C(D,oj)a\(Q,-u)) + 5&i 

The corresponding coefficients read: 
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(12) 


p 2V(t) ( lu H 6 I + l e l 7 E c) 


(13) 



where e = f2'/A defines the strength of 4WM, and: 



[/H = ^/[io;+| e p(ic;- 7go )] 2 +4| e p|n|2 ) (14) 

V(Lo)=Lu(uj + i~f sc )-\n\ 2 . (15) 



The terms S&s and <5d/ represent the field fluctuations 
corresponding to the Langevin noise operators. These 
can be neglected in the weak-field approximation since 
their normal ordered contributions are proportional to 
(<7 ee ) and (a ss ) respectively, which are second order in 
the signal field. However since these terms contribute to 
the noise, they will be considered in Sect. IV. 

Another important quantity for EIT based QM is the 
matter excitation, since the light field is mapped onto it 
during the storage process. This excitation is described 
by the spin operator a gs which can be found in terms of 
Eqs.(9,10): 



,(D,u) 



D,*g 
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as(D,uj) 
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(16) 



It has contributions from both the signal and idler field, 
but the idler term is proportional to the small parameter 
7 ge /A, and therefore can largely be ignored compared to 
the signal part. We can therefore concentrate on solving 
for the propagation of the signal field and our results will 
still be applicable to EIT quantum memory. Since we 
only consider propagation we ignore however any losses 
due to pulse leakage during the storing process. 

The coefficients A(D, to), C(D, w) describe the spectral 
transmission for the input signal and idler fields and the 
coefficient B(D,uj) describes the spectral coupling be- 
tween the fields. Since we are mainly interested in quan- 
tum memory applications and it has already been shown 
that an input idler field is not stored [18], we will assume 
no input idler field, implying we can ignore C(D,uS). In 
this case, the semi-classical solution for the signal and 
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FIG. 2: Plot of the transmission coefficient as a function 
of frequency for EIT with 4WM (solid) and without 4WM 
(dashed) . The parameters were chosen to emphasize the spec- 
tral behavior. 



idler fields are just: 

a s (D,to) 
ai(D,u) 



A(D,Lu)a s (0,u), (17) 

-B*(Z?, W )aJ(0, W ). (18) 



While the expressions for A and B are complicated, in 
the limit of £ < 1 which we have already assumed in 
order to derive the reduced Hamiltonian and D > 1 as is 
necessary for QM, both A and B are well approximated 
by Gaussians: 



A(D,u) = A (D)e iTs(D)u ^, (19) 

B(D,u) = B a {D)e Tl(D)ul ~~^i . (20) 



where the rg, tj are the group delay times for the signal 
and generated idler field, respectively; while Aujs and 
Aujj arc the frequency widths. It is clear that amplitudes 



Aq and Bq are the steady state solutions for the fields: 
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(21) 
(22) 



We introduce the ratio of the control field Rabi frequen- 
cies r\ = \Q!\/\il\, which only differs from unity when the 
two transitions have different dipole moments. It makes 
sense to introduce a parameter to keep track of the effec- 
tive 4WM strength, 



x = Drj 



A 



(23) 



Then for field propagation, Eqs.(21,22) defines two dis- 
tinct regimes of x. For large 4WM strength x > 1, both 
the signal and idler field experience exponential growth 
and except for a phase factor are essentially the same 
amplitude, as shown in Fig. 3. 
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FIG. 3: Log plot of the signal (solid) and idler (dashed) ampli- 
tudes as a function of effective 4WM optical depth. Assuming 
no initial idler field and normalized to the amplitude of the 
initial signal field. Note both become equal and grow ex- 
ponentially for high optical depths. The inset is a linear plot 
showing the the low optical depth behavior of the amplitudes. 



not to have a time delay. In the low optical depth case: 
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At large optical depths the 4WM process is generating 
many more photons than in the initial pulse, and for 
every new signal field photon there is a corresponding 
idler photon generated. While for smaller optical depth, 
I<1, there is only weak gain for the idler and signal 
field, i.e. we can treat this perturbatively: 



An = l 



X 



Bo = l im x - 



(26) 

(27) 



The idler field grows faster than the signal, but since it 
starts from vacuum, it remains much weaker than the 
signal field. Since 4WM introduces gain on the signal 
field, which also leads to stronger matter excitations, 
when present it will always increase the classical storage 
and retrieval efficiencies for an EIT memory. Therefore, 
in experiments that see a loss of classical efficiencies at 
higher optical depths such as [16], the loss should not be 
attributed to 4WM, but rather to other processes that 
grow with optical depth such as increased dephasing or 
depletion of the control field. The case of particular in- 
terest for QM will be for small 4WM strengths x < 1, 
since we will show in Sect. IV that exponential growth 
of the signal field is accompanied by an equally strong 
growth in noise. 

The group delay time for the fields is given by ts and 
r/, where a field traveling at the speed of light is taken 



such that ts is essentially the standard EIT delay time. 
The delay time for the generated idler 77, is approxi- 
mately half of that for the signal field. This can be under- 
stood as a consequence of the idler field being generated 
from the signal field. While the idler field is essentially 
moving at the speed of light and therefore is not delayed, 
it is also constantly being generated by the slow signal 
field. The total delay time is then the average of de- 
lay for idler photons generated near the beginning of the 
medium, and the idler photons generated at the end of 
the medium after the slower signal field has traversed the 
medium length. A similar effect is seen at high optical 
depth: 



TS — Tl 



2101' 



(30) 



where there is now a locking of the velocities for both 
fields to the average. Both fields are growing exponen- 
tially and travel together, the faster idler field is gener- 
ating a new slower signal field which leads to less signal 
delay, while the slower signal field generates the idler field 
leading to a longer idler delay than the low optical depth 
case. 

The spectral behavior of the fields described by the 
frequency widths Aws and Au>j, and we can again dis- 
tinguish between two different regimes. For low optical 
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FIG. 4: Plot of the signal delay time, normalized to the stan- 
dard EIT delay time .D7 ge /|fi| 2 (solid) and the signal trans- 
mission frequency width, normalized to the standard EIT 
transmission window |n| 2 /(7 gc \/i5) (dashed), as a function 
of the effective 4WM optical depth. 



depth, 




(31) 
(32) 



Where Atx>5 defines the usual EIT transmission window 
with a small correction, and Aw/ defines a larger window 
of transparency since the idler field does not experience 
absorption. At high optical depth: 



Aljs — Awj 



mi 



"/g<: 



8jgeV 



(33) 



The fields arc now propagating together with a similar 
transmission window that is narrower than the original 
EIT transmission window by a factor of y/&yg^/A. This 
narrowing is due to preferential gain of the signal near 
resonance where it does not experience absorption, rather 



than at frequencies near the edge of the EIT transmis- 
sion window that see some absorption. We have shown 
that for low optical depth D < |A|/7 gc , the signal field 
propagates similar to normal EIT, with a small gain due 
to 4WM from the newly generated idler field which re- 
mains weak compared to the signal field and propagates 
through the medium as if it was transparent. The prop- 
agation is dramatically different at high optical depths 
D > |A|/7 ge , where exponential growth of the signal and 
idler field couple the fields such that they have equal am- 
plitudes, experience less group delay, and have a narrower 
transmission window. 



IV. SIGNAL INTENSITY AND ADDITIVE 
PHOTON NOISE 

As an indicative measure for the effect of 4WM on light 
storage in an EIT medium we now consider the number 
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FIG. 5: Plot of the signal field at the beginning of the medium 
(dot-dashed), after an effective optical depth of i = .75 
(solid), and after x = 3 (dashed). Normalized to the initial 
pulse amplitude. 



of signal photons at the end of the medium at some time 
r using Eqs.(9,10): 



(al(D,T)a s (D,T)) = doj' du e'^' ~ u)t (A*{D,b/)^ 8 (p,u')A(D t u)& s (p i u)) 

+ {B*{D,uj')a I {D,Lo')B{D,uj)a\{D,uj)) 
= ffdu'duje-^'-^ {A*(D,u J , )al(0,uj')A(D,uj)a s (0,u J )} + f duj\B(D,Lj)\ 2 , 

I 



(34) 



(35) 



where the first part corresponds to the semi-classical solu- [ar(a/), a}(w)] = 5(u) — to'), i.e. this part is a pure quan- 
tion, which one would obtain by treating the fields clas- turn mechanical effect which does not appear in the semi- 
sically with no input idler field. For the second part classical solution. Since the value of the second part is 
we used the commutator relation for the field operators equal for all time r it describes the generation rate of the 

incoherent signal photons. This contribution exists even 
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FIG. 6: Plot of the number of noise photons produced as a 
function of effective 4WM optical depth. Notice it becomes 
larger than 1 slightly before when D > A/(27 go ). 



when there is no field input at all, and therefore is impor- 
tant for few photon input fields consequently, we refer to 
this as the vacuum noise contribution. It does not grow 
with signal field strength, so is much less important for 
fields with large photon number. With this we are able 
to estimate the number of noise photos by multiplying 
the generation rate with the propagation time of the sig- 
nal field. As noted in Sect. Ill, at high optical depth 
x > 1, Aq and Bq are equal, which implies that for a 
single photon input, the vacuum noise contribution will 
be as strong as the output of the signal field. Therefore 
in the regime of x > 1, a quantum memory is impossible. 
And as one can see from the graph in Fig. 6 we get addi- 
tional noise photons even for 4WM strengths not much 
larger than .5. 



V. FIDELITY OF PROPAGATION 

For a wave propagating through our medium, we can 
calculate the fidelity by finding the overlap of the wave- 
function well before the medium and the wave-function 
well after leaving the medium. 
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To describe propagation we consider a quantization 
length Lq much longer than the medium, which we sub- 
divide into 2M + 1 spatial intervals, each with center 
point: 



nL Q 
2M + V 



(37) 



with n running from — M to M. To describe a field pulse 
we also must divide frequency into 2M+1 discrete modes, 
centered at the resonance frequencies. We want to have 
two creation operators, one that creates a signal field 
photon of frequency ui n given by c„ , and one that creates 



an idler photon given by 6„, such that: 
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where we are labeling the vacuum for the signal and idler 

field as: 
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Then the field annihilation operator for the signal or idler 
field at a particular zi arc described by: 
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With the medium centered inside the quantization 
length, it begins at some z q and ends at z q - . Since when 
not in the medium the field is freely propagating, all 
we need to know is how it changes within the medium. 
Eq.(9) gives an input-output relation for each frequency 



J*-nCq,n 
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(44) 



which we can use to find the normal ordered characteris- 
tic function \n for a particular input state. The charac- 
teristic function in turn is the complex Fourier-transform 
of the multi-mode Glauber P-representation, from which 
we can easily determine the overlap in Eq.(36) defining 
the fidelity. 

We are particularly interested in the input of a single 
photon state in the signal channel, which we represent 
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where fi is a distribution function that represents the 
photons frequency envelope, and is normalized to pro- 
duce a single photon by: 
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After some straight-forward algebra we can then write 
the fidelity for a single photon input state as an integral: 
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where |({1}|{/3})| 2 is the overlap of the single photon 
state with the multimode coherent state: 

i({i}iw)i 2 =(II e " l/3 " |2 )(E/;/^^)' ( 48 ) 

n j,k 



and xn is the normally ordered characteristic function 
at the output of the medium which can be expressed 
in terms of the solutions of the Hcisenberg equations of 
motion for the field operators: 
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Now using Eq.(44) we get an expectation value that we 
can solve, in particular we can break the problem into 
finding the contribution from the signal field and the vac- 
uum contribution from the idler: 
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Eq.(53) has two parts, a product multiplied by a sum. 
The product can be interpreted as the vacuum contribu- 
tion, it is the same as would be calculated for a vacuum 
state. It always converges and depends on the spectral 
width of the coupling coefficient. In the continuous limit 
it can be explicitly calculated: 
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Performing the integrals while being careful with the 
sums gives the single photon fidelity as: 
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For large optical depths |i?o| 2 ^S> 1, this term dominates 
the single photon fidelity quickly dropping it to zero, with 
a rate that is at least exponential in optical depth. The 
second part is given by the sums in Eq.(53), and is due 
to the gain on the signal field; in the continuum limit 
M — >• oo can be converted back into integrals: 



E 



I f |2 A ' 
> Jl > V2AI+1 

(l + MIL) 



+ °°Mm? A ^ ] 



(i + I^MI 2 ) 



E 



E 



\l 



l + |Bo| 2 

12 l^.l 2 
1 V2AI+1 

Mil 



\ 



1 



i , A ^o o [go] 2 A "o 



exp 



V 1 T 2M+1- 

\f 12 l-Bjl 2 
l-'jl V2A/+1 



+oo 



dw|/(w)| 



l+|Bo| 2 Ac^ 

|A(c)| 2 



t^Acj^ 



IA)| 2 



1 



|J3ol 2 A^g 
l + \Bo\ 2 A^f 



(1 + |B( W )| 2 ) 1 + |Bo| 2 \ i + 2^4 - 2'?°'' £4 



, 2 ' 



l^l 2 



-^ Cl+ -L2ii_) 

3 \*-~ 2M+l> 



+oo 



<M/Mf 



|B(«) 



(l + |BM| 2 )-l + |Bo| 2 
I 



\ 



, l-|gp| 2 A^g 
i+ l + |Bn| 2 A^J 



(55) 

(56) 
(57) 



Where we have assumed that all the functions are Gaus- sians given by Eqs.(19, 20) and our initial distribution 
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FIG. 7: Plot of fidelity for the single photon case as a 
function of the effective 4WM optical depth. In the limit 
Ao/Aws, Acj /Au)i -> 0. 



function is a Gaussian such that: 
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(58) 



where Aloq is the frequency width of the incoming pulse. 
Now lets consider the ideal case where the frequency 
width of the incoming pulse fits well inside the transmis- 
sion window and define this ratio as y: 



Auj s 



(59) 



where we have chosen this limit such that the exponential 
term in Eq.(55) can be dropped. Then consider two dif- 
ferent limits for the 4WM strength, for x > 1 the fidelity 
exponentially decreases, while for small x and keeping 
only the first order terms we can express the fidelity as: 
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This shows us that the single photon fidelity will always 
be worse due to 4WM, although this is the expected re- 
sult since without 4WM, we have standard EIT, which for 
narrow pulse spectrum has fidelity close to 1. Therefore 
when 4WM is unavoidable, it would be best to imple- 
ment quantum memory in the regime D < A/j gc , which 
can be accomplished by choosing a medium with small 
7g C , such as picking a low temperature implementation 
like a cold Rubidium gas trapped in a magneto-optical 
trap, rather than a warm gas cell. 



VI. FIDELITY OF PROPAGATION WITH 
LOSSES 

As we have seen in the previous section, the presence 
of gain due to 4WM leads to a fast reduction of fidelity in 
otherwise loss-less propagation. We now analyze if 4WM 
could be beneficial when there are some linear losses due 
to scattering in the medium. We consider in particular 



the case where the 4WM gain exactly compensates linear 
loss, for which we will compare the EIT and 4WM single 
photon fidelities. 

Our analysis in Sect.V also applies to EIT in the limit 
of Aq — > 1 and Bq — > 0. In order to model linear losses 
in EIT with spatial loss coefficient A we just need to take 
the expression for the single photon fidelity and replace 
the coefficients with Aq = exp(— AD/2) and Bq = 0. In 
that limit the only integral needed is much simpler with: 
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which in the limit of Aujo/Aujs -C y/2/D, collapses the 
fidelity of EIT with losses to the expected result of: 
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(62) 



The same loss can be added to the 4WM fidelity by taking 
A Q -> A exp(-AD/2) and B Q -> B Q cxp(-\D/2). 

We pose the question of when does the 4WM sin- 
gle photon fidelity surpass that of the the EIT, assum- 
ing both systems experience linear loss. In the limit of 
Au>o/Au>s —> 0, we can approximate the fidelity as: 
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(63) 



Now we can take the steady state solution |^4o| 2 = 
qcosh 2 (x) and |B | 2 = <7sinh 2 (x), with the EIT fidelity 
given by q = exp(— XD). In this case it is fairly simple 
to calculate when 4WM can improve over EIT, it is pos- 
sible when q < .5, i.e. the EIT fidelity is already less 
than .5, as illustrated in Fig. 8. So while 4WM can be 
an improvement, it only helps in cases where the fidelity 
is already too bad to use as a quantum memory. 



VII. NOISE DUE TO FINITE EXCITED STATE 
POPULATION 

There will also be noise contributions due to sponta- 
neous emission, which we previously neglected since con- 
tributions from the Langevin noise operators are second 
order in the signal or idler fields. In the following we will 
calculate these noise contributions. The noise operators 
introduced in Eqs.(9,10) read explicitly: 
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FIG. 8: Plot of fidelity of four wave mixing (solid) and stan- 
dard EIT (dashed) fidelities as a function of effective optical 
depth in the presence of linear losses, under the assumption 
that the signal field is always well within the transmission win- 
dow. The losses are taken to match the 4WM gain present at 
the effective optical depth of 2, and with "/ ga g/A — 1/30. 



J 



where A, B, and C are given by Eqs. (11-13), and the new 
F operators are defined in terms of the Langcvin noise 
operators: 
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The rate of generation for extra amplitude noise due to 
spontaneous emission, is given by the expectation value 
(6a s Sas) for which to the first non-zero order in il'/A 
is: 
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While the form of the F are unknown, their expectations can be calculated using the fluctuation dissipation theorem 
[20], but we need to consider the full set of equations from the reduced Hamiltonian given in Eq.(l), not just the a ge 
and tJg S equations Eqs. (5,6). The full equations, neglecting detunings, letting 7 gs = 0, and dropping terms that are 
third power in the fields are: 
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Notice that the equations are no longer linear, but we 
expect the noise term to be second order in the fields. 
From Eqs. (70-73) we can find the expectation values for 
the Langevin noise operators: 



{£ ge^ gel 
i^gs^gs) 
\-P ge^ gsl 
X* 1 gs^gel 



(27gc - r es - r eg )(<r ee ), (74) 

resaee), (75) 

(7gc-7sc)(o-L), (76) 

(7gc-7 S c)(o- se ). (77) 



Now with these expectation values we can calculate 



(F'sFs): 
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The expectation values of the atomic operators can 
be found in terms of the field expectations by solving 
Eqs. (70-73), in the limit |7 so - 7 gc | < 7 SC the (a se ) 
can be neglected. Then only the excited state popula- 
tion is important. We can further simplify by assuming 
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w ^C 7go,7so, since this term is multiplied by the A(w) 
which vanishes for frequencies outside the EIT width, 
then the average excited state population can be found 
by solving the semi-classical form of Eqs. (70-73): 

_2g 2 lgc Lu 2 a s a s 2g 2 lgo u J (\n'n\/A)^/\asa I \ 



r eg \V(cj)\ 2 r eg \V(u)\ 2 
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where as and a] are the semi-classical field solutions 
given by Eqs. (17, 18). Replacing the initial field expec- 
tation value with the distribution given by Eq.(58) and 
multiplying by the delay time of the signal field tq yields 
the number of noise photons affected by dephasing due 
to spontaneous emission: 
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The integral over £ is straight forward when we use the Gaussian approximation developed in Sect. Ill for A(£, u>) 
and B(£,uj). Approximating |Ao(£)| = cosh( 7gc 7y£/A) and |£?o(£)l = smn (7gc^£/^)i and performing the £ integral 
leaves us with: 
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Now consider that the spectral window of A(D, to) will ensure that the main contribution of the integral comes for 
small u>, so we can take w < |f2| 2 / 7gc . This will allow us to approximate by replacing 7go w 2 + |n| 4 with |fi| 4 and 
perform the uj integral: 
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This can be simplified by taking r es = r eg , and by notic- 
ing that the square of the group index n 2 = g 4 A 2 /|f2| 2 
is approximately equal to l/( 7gG T£>). Assuming first that 
we are in the small optical depth regime, where x < 1, 
we can simplify further: 
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We can also consider the case when the optical depth is 
large x S> 1, then we have: 
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which like the vacuum noise is exponentially growing. It 
is important to note that the noise induced by sponta- 
neous emission is proportional to the initial signal field 



strength, and as shown in Fig. 9 for a small number of 
initial photons its contribution will be much weaker than 
the vacuum noise contribution given in Sect. IV. There- 
fore our dropping of this effect in the fidelity calculation 
of Sect. V is justified. Of course for a classical field with 
a large number of photons, dephasing due to spontaneous 
emission will be the dominate contribution to the noise. 



VIII. CONCLUSION 

We developed a model for a field propagating inside 
an EIT medium that has 4WM. We found that there 
are limits on the use of 4WM EIT as a single photon 
quantum memory based just on the propagation fidelity 
There arc two sources of noise, first wc identified a 4WM 
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FIG. 9: Number of photons affected by spontaneous emis- 
sion showing at what effective 4WM optical depth it becomes 
larger than 1, assuming a single photon input. For compari- 
son the vertical dashed line shows the effective 4WM optical 
depth when the vacuum noise becomes 1. 



process that creates new signal photons from the vacuum. 
The second source comes from the finite population that 
4WM adds to the excited state, leading to dephasing of 
the dipoles due to spontaneous emission. Together both 
sources of noise become exponentially large for optical 
depths D > A|£7'|/(7g C |S!|). This gives a natural limit 
on how high optical depth can be in EIT based quantum 
memories when 4WM is present. 

We further show that even in the best case scenario 
where the gain from 4WM compensates some natural 
losses in the system, for example due to scattering, the 
propagation fidelity of 4WM EIT is still worse than that 



for standard EIT unless the EIT fidelity is below one 
half. Therefore for an EIT quantum memory, it is always 
preferential to avoid four-wave mixing. This can be ac- 
complished by either choosing proper polarizations such 
that the control field can not couple to the signal tran- 
sition, or by keeping the optical depth below the limit 
D < A/j gc . Since high optical depth is needed for an 
efficient quantum memory, this limit for D is easier to 
accomplish in systems with smaller 7 ge , such as cold gas 
based memories rather than those in hot gas where the 
ratio of A/7 gc is much smaller. 

Our model is limited by only considering propagation 
of the field through the medium, this ignores the fact 
that it is actually the collective spin excitation that gets 
stored in an EIT memory. Our approach also neglects 
considering the limits imposed by needing the field to be 
wholly within the EIT medium at the point where the 
field would be turned off for EIT storage. While con- 
sidering these effects would not improve the limit 4WM 
imposes on the optical depth, it is possible that when 
considering the fidelity of the spin wave and the limits 
on the pulse length, that there are scenarios for lower 
optical depths where 4WM could be useful. We plan to 
further investigate the effect of 4WM on quantum mem- 
ory, in particular by finding the effect of 4WM on the 
collective spin state. 
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